Glass Transition in Confined Geometry 
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Extending mode-coupling theory, we elaborate a microscopic theory for the glass transition of liquids confined 
between two parallel flat hard walls. The theory contains the standard MCT equations in bulk and in two 
dimensions as limiting cases and requires as input solely the equilibrium density profile and the structure factors 
of the fluid in confinement. We evaluate the phase diagram as a function of the distance of the plates for the case 
of a hard sphere fluid and obtain an oscillatory behavior of the glass transtion line as a result of the structural 
changes related to layering. 
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Dense liquids display an intriguing complex dynam- 
ical behavior upon approaching the glass transition, as 
manifested for example in the drastic slowing down of 
transport, the appearance of stretched structural relax- 
ation, or the emergence of power laws at mesoscopic 
time scales 0]. For bulk systems, a scenario of the 
evolution of slow complex dynamics has been presented 
by the mode coupling theory (MCT) of the glass tran- 
sition, many aspects of which have been confirmed in 
the last two decades by experiments and computer simu- 
lations |Hl|2l. Yet, it is not obvious how the approxima- 
tions capture the collective rearranging of the local cages 
which is expected to occur via increasing cooperativity. 
Therefore significant research effort has been performed 
recently to confine the liquid to small pores, films, or 
tubes [3H5), since these experiments may hold the key to 
directly unravel the essence of the glass transition. 

Computer simulations have shown that walls, in gen- 
eral, induce already significant changes in the static 
structure and, in particular, smooth walls lead to lay- 
ering [6]. Second, it has been found that the dynami- 
cal features and transport properties quantified, e.g. by 
the diffusivity Q, are strongly influenced by confine- 
ment. Here, the liquid-wall interactions play an im- 
portant role: while Lennard- Jones interactions lead to 
higher glass transition temperatures by increasing the 
confinement (El|9), the opposite happens for purely re- 
pulsive walls iflOl . Interestingly, the drastic dynamical 
changes persist even in a model for an artificial pore, 
where the static structure remains identical to the bulk 
system (8][9). 

Experimental results from confocal microscopy ap- 
plied to colloidal suspensions between two quasi-parallel 
walls reveal a decrease of the critical packing frac- 
tion ifTTI due to confinement and a smaller mean-square 
displacement (MSD) parallel and close to the walls Ifl2ll . 
The MSD perpendicular to the walls exhibits an oscilla- 
tory dependence on the distance to the plates, in contrast 
to the MSD parallel to the walls IfTTI . A similar oscil- 
latory behavior of the diffusivities of a mono-disperse 



system of hard spheres with packing fraction cp = 0.40 
between two parallel and flat walls has recently been 
found from a computer simulation, however, for a di- 
rection perpendicular and parallel to the wall [7]. The 
glass transition of a single layer of binary colloidal mix- 
tures [13] displays essentially similar features as in bulk, 
which suggests that MCT is applicable also for the two- 
dimensional case 1 1414161 . 

In this Letter, we extend MCT to inhomogeneous liq- 
uids confined between two parallel flat hard walls with- 
out surface roughness. The derivation employs symme- 
try adapted eigenfunctions and a corresponding splitting 
of the current densities. Then the theory will continu- 
ously interpolate between the glass transition in two di- 
mensions and in bulk and relies solely on the equilibrium 
structure of the fluid in confinement. As an example, we 
consider a hard sphere fluid and determine the glass tran- 
sition line for various distances of the plates. In particu- 
lar, we obtain the oscillatory behavior of the glass tran- 
sition line as a result of the structural changes connected 
to the layering. 

Let us mention that MCT has already been applied to 
liquids in random porous media modeled by a quenched 
matrix of particles [17]. This kind of MCT is quite sim- 
ilar to MCT for the Lorentz model where a particle dif- 
fuses through randomly distributed obstacles |[T8l and 
has no relationship to what we aim at here. Another in- 
teresting approach is the extension of MCT to liquids in 
external fields [19] to calculate the linear response of the 
intermediate scattering function to a weak slowly vary- 
ing perturbation in order to deduce a divergent dynami- 
cal length. Here in contrast, we consider strong confin- 
ing potentials with variations even on microscopic length 
scales. As far as we know, MCT has never been extended 
to such situations, particularly to a slit, cylindrical ge- 
ometries, etc. where layering occurs. 

We consider a liquid of /V particles confined between 
two parallel flat hard walls which restrict the centers of 
particles to the slab between ±L/2. The microscopic 
density p(x,t) can be decomposed into suitable Fourier 
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modes 

N 
n=l 

where f n {t) - (x n (t),y n (t)) and z n {t) denote the position 
of the «-th particle parallel and perpendicular to the wall, 
respectively. The discrete values = 2nfi/L, /j. e Z ac- 
count for the confined geometry and q = (q x ,q y ) refers to 
the wave vector parallel to the walls. 

The walls induce a non-trivial density profile perpen- 
dicular to the confinement n(z) = (p(x,t)), that can be 
represented as a Fourier series with Fourier coefficients 

rL/2 

n M = J n(z)exp[iQ fl z]dz. For later purpose we also 
introduce the local specific volume v(z) := l/«(z); its 
Fourier coefficients are related to by a convolution 

The quantity of basic interest are the time-dependent 
correlation functions of density modes. Here, the inter- 
mediate scattering function is generalized to the infinite 
dimensional matrix S(q,t) = (S MV (q,t)) with elements 

Sfiv(qJ) = ■jz(6p tl (ef,tY8p v (if,0)), (2) 

where 6p M (q,t) = p M (q,t) - (p M (q,t)) denotes the density 
fluctuation. Mirror reflection and time inversion symme- 
try (for Newtonian dynamics) implies n~ M = n M , V-^ - v^, 
and S -ft- v (q, t) = S ^(q, t) = S v ^{q, t) = S MV (q, t)* , in par- 
ticular S fiv{q,t) is real symmetric. 
The continuity equation relates the densities p M (q,t) to 

a corresponding current j^iqj). The crucial step, in con- 
trast to bulk liquids, is now to decompose into a compo- 
nent parallel and perpendicular to the wall, a € {||,-L}, 

N 

= 2 b a (l f n (t), z n (t)) expUQ^n®] e^-W . (3) 

n=\ 

Here, we abbreviated q = q/q, and introduced the selec- 
tor b a {x,z) = x6 a> \\ +z5 a ,x- 

Choosing 6p M (q,t), j^(q,t) and j^(q,t) as distin- 
guished variables, the Zwanzig-Mori projection operator 
formalism E01I2TH leads to 

S(q,t) + f dt'K(q,t-t')S~ l (q)S(q,t') = 0, (4) 
Jo 

where S(q) = S(q, t - 0) is the static correlation function. 
The memory kernel K(q, t) is decomposed according to 

<*v8=±,|| 

and its components 'Kiq, t) - ( ( K^(q, t)) satisfy 

7C(q, t)+ f dt'J(q)M(q, t - f)7C(q, t') = 0. (6) 
o 



The matrix of the static current density correlators 

ST(q) - (3/lv(q)) inherits an explicit dependence on the 
mode index via the average density profile 

3%{q) = (Wffy®) = ^ (7) 

where m denotes the mass of the particles. 

Following MCT for bulk liquids [1] the memory ker- 
nel matrix M(q, t) = (M^y(q, t)) can be approximated as 
afunctional of the density correlators 

M$(q,t)*9$[S(k,t);q\= (8) 
Yj ^rt"iP2;vviv 2 (9> < 7i> < 72)*S' nm(<ll,t)S ^ 2 v 2 072,0> 

with the vertices 

x^v rt [fc tt (f 91, Q K ^ 2 )c M , ^ (q x ) + (1 <-> 2)] X 

k,k' 

x [^•<?i,G^-v 2 )c yi y_ V2 (<?i) + (l «-» 2)]v^_ v . (9) 

Here, the direct correlation functions c(q) - (c^ v (q)) are 
related to the static correlators S ^(q) by the inhomoge- 
neous Ornstein-Zernike equation l22l 

S- I (q)=j^[v-C(q)], (10) 

and (y) MV = v M . v . 

The MCT approximation for the memory kernel leads 
to a closed set of Eqs. ((4} - ((9]) for the density corre- 
lators SfjyiqJ). Note that this set is valid for any one- 
component liquid between two parallel flat hard walls 
with arbitrary particle-particle and arbitrary particle- 
wall interactions. These interactions only enter into the 
equations via the static quantities n^,v v and S ^ v {q). Let 
us emphasize that in contrast to bulk systems also the 
Fourier coefficients v M of the inverse density profile enter 
the vertices. One can prove [23 ] that this set reduces to 
the well-known MCT equations for a bulk liquid (2 for 
L — » oo and to those derived for a two-dimensional liq- 
uid, L -> 0, (HI- Consequently the MCT equations Q 
- (|9]> interpolate between a two- and a three-dimensional 
liquid, as required. Since the glassy dynamics is iden- 
tical for both Newtonian and Brownian motion [H, our 
results apply also to colloidal suspensions. 

To locate the glass transition point one has to solve 
the self-consistent set of equations for the non-ergodicity 
parameters F^ v (q) := lim^ooS MV (q,t): 

F{q) =[S-\q)+S- l (q)K(q)S- l (q)T l , (11) 
(K(q)) lty = Y b a (q,Q M )(T- l [F(ky,q]) l fy(q,Q v ). 
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FIG. 1: Static correlators Soo(q) for <p = 0.42 for different distances 
H of the plates. The first sharp diffraction peak varies nonmonoton- 
ically; lowest for H = 2.0<x, highest for H = 2.5a and intermediate 
for H = 3. Oct. Inset: The corresponding off-diagonal structure factor 
S()\(q) is negative and becomes smaller upon increasing H . 



In the following we determine the glass transition line 
for hard spheres of diameter cr confined between two 
hard walls of distance H = L + cr. For the numerical so- 
lution of Eq. ( fTT) ) one has to truncate the Fourier indices 
\/u\ < M and q has to be discretized |[23l . Due to the large 
increase in complexity the maximum possible value has 
been M = 3. In the two-dimensional limit H/cr — » 1 only 
the lowest mode contribute, S fi V (q,t) — > Soo(<7,f)<^Vc>vO- 
Already for H/cr > 2 higher order modes > have to 
be taken into account, particularly for H/cr » 1 . There- 
fore we have also performed a subsequent diagonaliza- 
tion approximation (DA), where only the diagonal ele- 
ments S m {q), c MM (q), F m (q) and T^{{F flli {k)};q\ have 
been taken as nonzero, allowing to handle M = 10. 

The static input quantities were obtained by first cal- 
culating the density profile n(z) within fundamental mea- 
sure theory EM . These results are used as input into the 
inhomogeneous Ornstein-Zernike equation which is then 
solved with the Percus-Yevick closure relation [25]. The 
depletion forces give rise to a strongly enhanced density 
profile in the vicinity of walls, which results in the char- 
acteristic layering. The principle peak in the static struc- 
ture factor S oo(^) at qcr « 2n is higher for the half inte- 
ger distance H/cr = 2.5 than for both integer distances 
H/cr = 2.0 and 3.0, see Fig. [TJ Consequently the static 
correlations appear reduced in the case that a few mono- 
layers just fit between the walls, and one anticipates this 
layering effect to have drastic influences on the glass 
transition. 

The maximum packing fraction <p := nncr 3 /6 (n = 
no/H), for which we were able to reliably calculate the 
static input quantities has been for <p = 0.42, so far, not 
large enough to induce a transition. Therefore we mimic 
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FIG. 2: Critical parameter Xc(H;tp = 0.42) for M = 3, with DA and 
without DA, and the dimensionless diffusivities Oy and D ± , respec- 
tively, taken from |7] for <p = 0.40. 



the increase of the static correlations by an additional 
multiplicative control parameter x m front of the MCT 
functional T^y ■ Let (f c (H) be the critical packing frac- 
tion for x — 1 ■ Then, if tp is below and close to ip c (H) 
the vertices vary linearly with <p, and we therefore ex- 
pect that the //-dependence of <p c (H) is qualitatively well 
described by the //-dependence of the critical parame- 
ter Xc(H;ip). Figure [2] depicts Xc(H; 0.42) obtained with 
and without DA for M = 3; both graphs strongly resem- 
ble, particularly in the region 2.0 < H/cr < 4.0 where Xc 
has larger variation in contrast to H/cr > 4.0. The larger 
absolute values of Xc without DA may originate from a 
partial cancelation in the vertices. Increasing the cut-off 
M in the DA, i.e. taken into account higher order modes, 
we found that Xc(H; 0.42) does not change qualitatively 
(not shown). 

As argued above, we expect for (p c (H) qualitatively a 
similar //-dependence as we have found for Xc(H;0A2). 
Now assume tp fixed below <p c (H). Then, the dis- 
tance (tp c (H) - tp) to the glass transition point is maxi- 
mum (minimum) for those H for which the Xc(H;<p) is 
maximal (minimal). MCT predicts for the diffusivity 
D MCl (H,<p) ~ [tp c (H) - (pyW where y(H) does not de- 
pend sensitively on the type of the liquid and is close 
to two. Therefore assuming only a weak dependence 
of y on H, MCT predicts a maximum diffusivity at 
(///o-)Mf x T w 2.0, « 3.1, « 3.8 and « 5. These values are 
surprisingly close to the corresponding values from the 
simulation, Fig. [2J particularly for (D ± ) sim . Note that 
the variation of the diffusivity with H is much less pro- 
nounced for H/cr > 4.0 than for H/cr < 4.0 which is re- 
produced by MCT, as well. 

In conclusion, we have generalized MCT to liquids 
confined between two parallel flat walls for arbitrary 
particle-particle and particle-wall interactions. The the- 
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ory differs in several important aspects from correspond- 
ing one in bulk and in the plane. First, the walls ex- 
change momentum with the fluid implying that the in- 
termediate scattering function is no longer diagonal in 
the wavenumber. Second, the theory requires a splitting 
of the currents in order to capture the limits H — » cr and 
H — » oo. Furthermore the theory requires explicitly the 
density profile induced by the confinement, in addition to 
the generalized structure factors, which renders the glass 
transition line sensitive to layering effects. We have ex- 
emplified the predictive power of the theory by repro- 
ducing the non-monotonic behavior of the diffusivity for 
hard spheres on the confinement. Particularly, its differ- 
ent behavior for integer and half-integer values of H/cr 
can be interpreted as follows. If H/cr equals an integer r 
then exactly r monolayers fit into the slit. If the mono- 
layers were perfectly flat, they could slide. In order to 
get a structural arrest tp has to be increased. If H/cr de- 
viates stronger from r the monolayers become rougher 
and penetrate into each other which reduces the parti- 
cle's mobility and hence the packing fraction for a struc- 
tural arrest. This interpretation is consistent with the H 
dependence of the principal peak of Soo(q), see Fig.[T] 

The //-dependence of D has also been related to that 
of static quantities such as the excess entropy s ex and the 
available volume fraction ~pq. Knowledge of the bulk dif- 
fusivity Db u ik then yields D(H, ...) ~ DbuikW//), - - •) for 
x - s sx ,~po provided D is not too small Boll . In MCT x 
corresponds to <p c {H) and, if y{H) is almost constant, a 
similar relation D MCT (//, . . .) « D^((p c (H), . . .) holds. 
Calculating s ex (H) and ~pq{H) from fundamental mea- 
sure theory, we find that s ex and J5o are not constant on 
the glass transition line. However, the H dependence of 
(p c is in phase with that of s ex and ~pq demonstrating the 
importance of the latter for the glass transition. 

Recently, it has been shown that bulk binary mixtures 
of hard spheres E71 l28ll or hard disks lPT5l display an 
intriguing dependence on composition. In confinement 
these mixing effects compete with the layering induced 
by the walls and should give rise to an even richer phase 
behavior. It is instructive to also compare the glass tran- 
sition line with the random close packing values. For 
hard disk mixtures in the bulk both lines exhibit a strik- 
ing similarity 031, and the random close packing of a bi- 
nary mixture of hard spheres in confinement displays an 
oscillatory dependence on the distance of the walls [29]]. 
Thus, it would be interesting to extend MCT also to mix- 
tures in confinement. In addition, the theory presented 
here constitutes a promising framework to rationalize the 
numerous glass-forming systems in simple confinements 
that have been investigated in the recent past Q. Fur- 
thermore, it is a challenge for the future to extend MCT 
also to pore models for which the static properties are un- 
effected whereas the dynamical features sensitively de- 
pend on the confinement El . 
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